Constraining dark energy with cross-correlated CMB and Large Scale Structure data 
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We investigate the possibility of constraining dark energy with the Integrated Sachs Wolfe effect 
recently detected by cross-correlating the WMAP maps with several Large Scale Structure surveys. 
In agreement with previous works, we found that, under the assumption of a flat universe, the ISW 
signal is a promising tool for constraining dark energy. Current available data put weak limits on 
a constant dark energy equation of state w. We also find no constraints on the dark energy sound 
speed c\. For quintessence-like dark energy (cj! = 1) we find w < —0.53, while tighter bounds are 
possible only if the dark energy is "clustered" {c 2 e = 0), in such a case —1.94 < w < —0.63 at la. 
Better measurements of the CMB-LSS correlation will be possible with the next generation of deep 
redshift surveys. This will provide independent constraints on the dark energy which are alternative 
to those usually inferred from CMB and SN-Ia data. 

PACS numbers: 98.70.vc,98.80.Es 



I. INTRODUCTION 

The recent Wilkinson Microwave Anisotropy Probe 
(WMAP) satellite measurements of the Cosmic Mi- 
crowave Background (CMB) anisotropy spectra provide 
an accurate determination of several cosmological param- 
eters 0. In combination with complementary results 
from galaxy surveys [3, t| these data strongly suggest 
that the present energy budget of the universe is domi- 
nated by an exotic form of matter which is also respon- 
sible for the present phase of acceleration as directly in- 
ferred from measurements of luminosity distance to Su- 
pernova type la 0. 

The presence of a cosmological constant term A in Ein- 
stein's equation of General Relativity (GR) is the sim- 
plest explanation for dark energy. Remarkably the A 
cold dark matter scenario (ACDM) is the minimal model 
to consistently account for all observations and therefore 
has emerged as the standard model of cosmology. 

The WMAP temperature anisotropy maps have been 
also cross-correlated with several surveys of Large Scale 
Structure (LSS) distribution and a positive correlation 
signal has been detected [U IS H This is indeed 
another success of the ACDM model since the existence 
of a positive correlation was already predicted nearly a 
decade ago by Crittenden and Turok as a test of flat A- 
cosmologies [l0|. 

Despite the simplicity of this concordance model, the 
nature of dark energy is far from being understood. In 
fact the existence of a small non-vanishing cosmological 
constant A arises questions to which present theoretical 
particle physics has found no consistent explanation. On 
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the other hand more exotic forms of matter such a scalar 
field cannot be a priori excluded. They have been pro- 
posed in several context as candidate for dark energy 
IH 13 13 El Eland found to be consistent with cur- 
rent observations Jig. Alternatively in a number of sce- 
narios it has been suggested that what appears as a dark 
energy component could be a manifestation of a modified 
theory of gravity or a consequence of spatial extra dimen- 
sions E3 ■ I n this paper we will consider only GR gravity 
with a dark energy fluid whose energy momentum tensor 
violates the strong energy condition. 

Within this framework the dark energy and the cosmo- 
logical constant differs for two main aspects: the latter 
behaves as a homogeneous fluid with a constant energy 
density, while the former is a non-homogeneous fluid with 
a time dependent energy density and pressure. A simple 
way of describing these models is by specifying the equa- 
tion of state w = pde/ Pde, where pde is the pressure 
and pde is the energy density. The cosmological constant 
corresponds to the specific constant value w = — 1, while 
a general dark energy fluid may have a time dependent 
equation of state w(t). General covariance requires A 
to be a homogeneous component, while dishomogencitics 
may occur in fluids with w ^ —1. The clustering proper- 
ties of different dark energy models can be parameterized 
by an effective sound speed defined as the ratio between 
the pressure to density perturbations in the rest frame of 
dark energy, (? e — Spde/Spde E3- Since we lack of a 
consistent theoretical formulation of dark energy we are 
left with constraining some phenomenological motivated 
form of w(t) and (? e . 

The cross-correlation between CMB and LSS offers a 
new complemen tary way of constraining these parame- 
ters 0, 120I I2H 1221 |23| . In fact the correlation between 
these data sets is consequence of the Integrated Sachs- 
Wolfe effect [24|, which has been shown to be a sensi- 
tive probe of the evolution and clustering of dark energy 
[25I |29 | . It is therefore timely to investigate whether cur- 
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rent measurements of the CMB-LSS correlation can al- 
ready provide some novel constraints on the dark energy 
properties. 

The paper is organized as follows : in Sec.[H]we intro- 
duce the ISW-correlation and study its dependence on 
dark energy. In Sec. II III we describe the data analysis. 
In Sec. lIVI we discuss the results and finally we draw our 
conclusions in Sec. 

II. THEORY 

In a flat dark energy dominated universe the gravi- 
tational potentials associated with the large scale struc- 
tures decay as consequence of the accelerated phase of 
expansion. CMB photons which cross these regions ac- 
quire a shift which generates temperature anisotropies. 
This is the so called Integrated Sachs- Wolfe (ISW) effect 
[24|. A natural consequence of this mechanism is that if 
there is a clump of matter, such as a cluster of galaxies 
in a given direction of the sky, we are likely to observe a 
spot in the corresponding direction of the CMB provided 
that the CMB photons have crossed that region during 
the accelerated epoch. We therefore expect to measure a 
positive angular correlation between CMB temperature 
anisotropy maps and surveys of the large scale structures. 

The ISW temperature fluctuation, Aisw: m the direc- 
tion 71 is given by: 

A ISW (ii) = -2 J e- T ^ — (i u z)dz, (1) 

where $ is the Newtonian gauge gravitational potential 
and e~ T< - z > is the visibility function to account for a pos- 
sible suppression due to early reionization. 

The density contrast corresponding to a clump of lumi- 
nous matter observed by a given survey in the direction 
72 is: 

dLSsdz) = b J (j>(z)6 m (i 2 ,z)dz, (2) 

with 5 m the matter density perturbation, b the galaxy 
bias and <j)(z) is the selection function of the survey. 
The 2-point angular cross-correlation is defined as 

C x (8) = (Aiswd^Lssh^)), (3) 

where the angular brackets denote the average over the 
ensemble and 6 = |t"i — 72 1. For computational purposes 
it is convenient to decompose C x (9) into the Legendre 
series such that, 

00 0/11 

C X (e)=J2^r^P l (cos(9), (4) 
1=1 

where Pi (cos 9) are the Legendre polynomials and C x is 
the cross-correlation power spectrum given by |22l |23| : 

C? = A*l-J d ^Alll sw {k)lt SS {k), (5) 



where is the primordial power spectrum. The inte- 
grand functions lf sw {k) and I^ ss (k) are defined respec- 
tively as: 

l[ SW (k) = -2 J e-^^Mkriz^dz (6) 

I^ ss (k) = b J ^SiizMkr^dz, (7) 

where $k and J™, are the Fourier components of the grav- 
itational potential and matter perturbation respectively, 
jl\kr{z)\ are the spherical Bessel functions and r{z) is the 
comoving distance at redshift z. 

Direct measurements of the cross-power spectrum C x 
are more robust for likelihood parameter estimation since 
these data would be less correlated than measurements 
of C x {9). However current observations provide only 
estimates of the angular cross-correlation function, for 
this reason we focus on C x (9). 

In order to better illustrate the model dependence of 
the ISW we compute the angular cross-correlation func- 
tion for different values of w and c 2 . We first compute the 
cross-power spectrum of a given model using a properly 
modified version of CMBFAST code [2^ , then we evalu- 
ate C x (9) using Eq. In this way the monopole and 
dipole contribution to the angular correlation function 
are subtracted by construction. We have implemented 
the dark energy perturbation equations as described in 
[2^| . In what follows we assume a flat universe with Hub- 
ble parameter h — 0.68, Qde = 0.7, baryon density 
ftbh 2 — 0.024, scalar spectral index n s — 1, no reion- 
ization. The amplitude of the primordial fluctuations 
A s = 0.86 (as defined in [29|) and the bias 6 = 1. We 
consider a Gaussian selection function peaked at z = 0.1 
with variance a<p = 0.15. 

In Figure we plot the angular cross-correlation for 
the case c~ = 1 and in Figure QJd for c 2 — 0. The different 
lines correspond to w = —0.8 (solid), w — —0.4 (long 
dash-dot) and w = —4 (short dash). 

We notice that at small angles the angular cross- 
correlation is characterized by a nearly constant plateau, 
while it rapidly vanishes at larger angles (8 > 10°). The 
overall amplitude of the signal up to angles 9 ~ few° 
is particularly sensitive to both w and c 2 .. This can be 
better seen in Figure QJ; and Figure where we plot the 
amplitude of cross-correlation at the plateau, C x (0.3°), 
as function of w in the case c 2 = 1 and c 2 — respec- 
tively. 

For Cg = 1, the amplitude has a maximum around 
w = — 1 and slowly decreases for decreasing values of w, 
while it rapidly falls to zero for w — ► (Fig-QJ;). This is 
because the dark energy contribution to the ISW effect 
is mainly due to the background expansion. In fact for 
models with w > — 1, as it) — >0 the dark energy driven ex- 
pansion is less accelerated and tends to the matter dom- 
inated behavior. Hence the variation of the gravitational 
potentials is smaller and consequently produces a neg- 
ligible amount of ISW as w — » 0. Similarly for models 
with w < —1, the dark energy affects the expansion later 
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FIG. 1: Angular cross-correlation function for = 1 (a) 
and Cg = (b), the different lines corresponds to w — — 0.8 
(solid), w — —0.4 (long dash-dot) and w = —4 (short dash). 
Amplitude of the angular cross-correlation at the plateau as 
function of w for cj, — 1 (c) and Ce = (d). 



than in models with w > — 1. This effectively extends 
the period of matter domination which leads to a lower 
ISW signal. 

Since a smaller ISW signal can be compensated by in- 
creasing the amount of dark energy density Qde, we ex- 
pect a precise degeneracy line in the £Ide — w plane. In 
particular lower negative values of w will be counterbal- 
anced by higher values of FIde- 

On the contrary for (? e = 0, the amplitude of the 
cross-correlation is a monotonic decreasing function of w 
(Fig. [TJl). In this case the decay of the gravitational po- 
tential is sensitive to the clustering of dark energy which 
is more effective as w decreases. Thus the amplitude of 
the ISW increases as w decreases jl^, • We therefore 
expect the degeneracy in the £Ide — w plane to be or- 
thogonal to the previous case. In fact increasing £Ide 
will compensate for larger values of w. 

This trend hold independently of the selection function 
as long as it is centered in a range of redshifts up z ~ 
0.7 — 0.8 for models with w > —1. However one might 
expect this to not be the general situation in the case of 
dark energy models with a time dependent equation of 
state |23J. 

In Figure |21 we plot the amplitude of the angular cross- 
correlation at the plateau as function of the redshift. We 
have used a Gaussian function centered at different red- 
shifts in steps of 0.1 and with a width = 0.05. The 
upper (lower) panel shows the case = 1 = 0), the 
different lines correspond to the same models as in Fig- 
ure^ It can be seen that for z > 0.2 the signal decreases 



FIG. 2: Amplitude of the angular cross-correlation at the 
plateau as function of the redshift. A Gaussian selection func- 
tion with width cr^ = 0.05 has been used. The different lines 
correspond to models as in Figure 



with the redshift in way that is strongly dependent on 
the dark energy parameters. Therefore redshift measure- 
ments of the cross-correlation are a potentially powerful 
tool to distinguish between different dark energy models. 

We may also notice that having used a narrower selec- 
tion function the amplitude of the signal is systematically 
smaller than that found in previous cases. For instance 
for w = -0.8 and c\ = 1 we have C(0.3°) = 0.06fj,K 
at z = 0.1 (FigJ2J upper panel), which is smaller than 
C(0.3°) = O.UfiK (see FigUJi) obtained using a wider 
selection function (er^ = 0.15). 

Hence a sharper selection function gives a smaller 
cross-correlation signal, eventually leading to larger un- 
certainties. On the other hand increasing the number 
of uncorrelated redshift bins would allow a better recon- 
struction of the redshift evolution of the cross-correlation. 

This suggests that there could be an optimal way of 
designing a large scale structure survey which maximizes 
the discrimination power of the cross-correlation inde- 
pendently of the model of dark energy. This can be 
achieved using Integrated Parameter Space Optimization 
(IPSO) techniques |2Jj. We leave this interesting possi- 
bility to further investigation and we refer to [22j for a 
detailed analysis of dark energy parameter forecast from 
future cross-correlation data. 



III. METHOD AND DATA 

We perform a likelihood analysis usin g th e collection 
of data presented in Gaztanaga et al. |33j, which has 
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the advantage of being publicly available and easy to im- 
plement. These data consist of measurements of the av- 
erage angular cross-correlation around 8 = 5° between 
WMAP temperature anisotropy maps and several LSS 
surveys. The angular average around 9 = 5° ensures 
that the signal is not contaminated by foregrounds such 
as the SZ or lensing effects which are relevant at smaller 
angles (8 < 1°). Taking the average does not wash out 
the sensitivity on the dark energy parameter since, as 
we have shown in Section [H] the amplitude of the cross- 
correlation around 8 = 5° is still strongly dependent on 
the value of w and c 2 e . 

Possible systematic contaminants, such as extinction 
effects, seem not to affect these data and for a more de- 
tailed discussion we refer to . 

The data span a range of redshift O.f < z < 1 and 
for each redshift bin the data include an estimate of the 
galaxy bias with 20% errors. These biases are inferred 
by comparing the galaxy-galaxy correlation function of 
each experiments with the expectation of LCDM best fit 
model to WMAP power spectra. 

The amplitude of the scalar density perturbation A s is 
an overall normalization factor which we can marginal- 
ize over, on the contrary prior knowledge of the bias is 
required. In fact a scale and/or redshift dependent bias 
can in principle mimic the redshift evolution of the cross- 
correlation predicted by different dark energy models. It 
is therefore necessary to have an independent estimate 
of b, for instance by combining weak lensing information 
|34j or measurements of the matter power spectrum with 
the bispectrum |3fij . 

One of the advantages of testing dark energy with the 
cross-correlation is that it is insensitive to other param- 
eters which limits common dark energy parameter ex- 
traction analyses involving CMB temperature and po- 
larization anisotropy spectra |3fij |. For instance the ISW 
correlation is not affected by a late reionization or by an 
extra background of relativistic particles which c hang e 
the CMB spectra through the early-ISW (see e.g. |37j). 
The ISW-correlation is also independent of the ampli- 
tude of tensor modes and depends uniquely on the scalar 
perturbations, since a primordial background of gravity 
waves is uncorrelated with present large scale structure 
distribution. 

There is little sensitivity to the scalar spectral index 
n s , while the dependence on the baryon density Q& can 
be non-negligible. In fact the presence of baryons in- 
hibits the growth of CDM fluctuations between matter- 
radiation equality and photon-baryon decoupling [38j | 
causing the matter power spectrum to be suppressed on 
scales k > k eq for increasing values of Qb (k eq is the 
scale which enters the horizon at the equality) . Over the 
range of scales which contribute to the ISW-correlation 
(k ~ O.Of) the sensitivity on f2& is still present. In 
order to limit the number of likelihood parameters we 
therefore assume a Gaussian prior on the baryon density 
fl b h 2 = 0.0216 ± 0.002 consistent with WMAP and Big- 
Bang Nucleosynthesis bounds and consider the following 



set parameters: the matter density f2 m , the Hubble con- 
stant h, the equation of state w and the sound speed of 
dark energy c 2 .. We assume a scale invariant primordial 
spectrum n s — 1 and fix the optical depth to WMAP 
best fit value r = 0.17 (again the ISW is not particularly 
affected by a change in those parameters). We marginal- 
ize over the normalization amplitude A s , although we 
found no difference assuming the WMAP value. In fact 
changing A s shifts the overall amplitude of the angular 
cross-correlation of the same amount over different red- 
shifts but it does not change the redshift dependence of 
the signal. 

Since the experimental data are corrected for the bias 
by comparing the measured galaxy-galaxy correlation 
function in each redshift bin to the WMAP best fit model, 
we rescale these biases to each of the dark energy model 
in our database as described in |3flj . 

We compute for each theoretical model the angular 
cross-correlation as described in Sec. |n] using the selec- 
tion function 

^(z)~z 2 exp[-(z/z0 1 - 5 ], (8) 

where z% is the median redshift of the i-th survey. Then 
following |33| . we compute the average cross-correlation 
in the i-th bin, Cf- , around 8 ~ 5°. 

The data points are an average over angles and are in- 
ferred from surveys whose selection functions may over- 
lap in redshift space, hence they are not independent 
measurements and indeed are affected by a certain de- 
gree of correlation. Since we have no access to the raw 
data we have no way of accounting for the first type of 
correlation, while using Eq. (JHJ) we can estimate the cor- 
relation between different redshift bins. We compute the 
correlation matrix p = {py}, where pij is the fraction 
of overlapping volume between the i-th and j-th sur- 
veys (i.e. the diagonal components are pa = 1). To 
be more conservative we have assumed that the different 
surveys cover the same fraction of sky, in general this 
is not the case and the fraction of overlapping volume 
can be smaller. We found that only two data points are 
highly correlated, since their selection functions overlap 
for about 70% in redshift space, while the correlation 
among the remaining data points are less than 20%. 

We compute a likelihood function £ defined as 



- 21og£ = x 2 = ]T(Cf - Cf)Mr. 1 (Cf - Cf), (9) 
y 

where Cf are the data and MZ = pij/(<Ji<Tj) is our 
estimate of the inverse of the covariance matrix, with <jj 
the measured uncertainty in the i-th bin. 

We also use the SN-Ia "Gold" data [3^| to derive com- 
plementary constraints on the dark energy parameters 
and compare with those derived from the ISW correla- 
tion. 
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FIG • 3: Two-dimensional marginalized likelihoods on Q.de — 
w. The yellow and red area correspond to 1 and 2a limits 
inferred from the ISW data for = 1. Solid and dash lines 
represent the 1 and 2er contours from the SN-Ia data. 



IV. RESULTS 

We now discuss the results of the likelihood analysis 
and show the marginalized constraints on the dark en- 
ergy parameters. As expected we found that the results 
depend on the dark energy clustering. For instance in 
Figure and 0] we plot the two-dimensional marginalized 
likelihoods in the Ude — w plane for c 2 = 1 and c 2 = 
respectively. The yellow and red contours correspond to 
1 and 2a limits. We also plot the marginalized likelihood 
inferred from the SN-Ia data. 
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FIG. 4: As in Figure Hwith cl = prior. 

It appears evident that the derived limits on w strongly 
depend on the prior value of . In particular varying c 2 e 
changes the direction of the degeneracy between CIde 
and w. This is a natural consequence of the discussion 
presented in Section [D] Namely, for c 2 = 1 decreasing 



w causes a lower ISW signal which can be compensated 
by larger values of SIde, while the opposite occurs for 
c 2 = 0. Therefore for c 2 = 1 the ISW data provide only 
a weak upper limit on the equation of state, w < —0.53 
at 68% confidence level 43] . On the contrary assuming 
c 2 = we obtain a tighter constraint —1.94 < w < —0.63 
at 2a. 

We can also notice that for c 2 = the combination 
of the ISW likelihood with the SN-Ia one has little ef- 
fect in determining w, since the Qde — w degeneracy of 
the luminosity distance lies in the same direction of the 
ISW cross-correlation. This is not the case for c 2 = 1, 
where the degeneracy lines are orthogonal. Therefore 
the combined likelihoods provide a more stringent bound, 
-1.72 < w < -0.53 at 95%. 

These limits are stable assuming a standard prior on 
the matter density. For instance for Cl m = 0.27 ±0.04 we 
obtain —1.51 < w < -0.72 for c 2 e — and -1.81 < w < 
-0.53 for (? e = 1 at 95%. 

The constraints we have inferred so far would be 
15 — 20% tighter if we had ignored the correlations 
amongst different redshift bins. In fact, as we have dis- 
cussed in the previous Section, the level of correlation 
of current data is still non-negligible. On the contrary 
the next generation of surveys will be characterized by 
more localized selection functions and provide uncorre- 
cted cross-correlation measurements. 

In Figure [S] we plot the 5 data points and three dif- 
ferent model predictions of the average angular cross- 
correlation normalized to the bias and the data. The solid 
line corresponds to the LCDM best fit model, we also plot 
a non-accelerating dark energy dominated model (long- 
dash line) and a phantom model (short-dash) which are 
disfavored by the data. 




mean redshift (z) 



FIG. 5: Dots with errorbars are the different measurements 
of C x /b, the solid line corresponds to the LCDM best fit 
model, the long-dash and short-dash lines show the case (w = 
—0.2, Q,de = 0.7) and (w = —2, Q DE = 0.6) respectively. 
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In Figure|B]we plot the marginalized 1 and 2a contours 
in the log 10 c\—w plane, as it can be seen the dark energy 
sound speed remains unconstrained. This is consistent 
with the results from the CMB data analysis by Weller 
and Lewis [3(J and only at 2a with Bean and Dore psj . 
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FIG. 6: Two-dimensional marginalized likelihood on log 10 c e — 
w from ISW-correlation data. The yellow and red contours 
correspond to 1 and 2a limits. 




FIG. 7: 1 and 2a limits from ISW-correlation (yellow and 
red contours) and SN-Ia luminosity distance (dash and solid 
lines) on too — »i . 

We have also extended our analysis to constrain a class 
of slowly varying dark energy models described by a pa- 
rameterization of the equation of state which linear in 
the scale factor: w(a) = wq + (1 — a)wi [ill ]. 

This part should be considered as a simple exer- 
cise since current cross-correlation data are not accurate 
enough to allow us to constrain more than two param- 
eters. We therefore assume an Q m = 0.3 prior and we 
also restrict our analysis to models satisfying the Weak 
Energy Condition (WEC) , wq > — 1 and < w\ < 1 such 



that we can consistently account for the dark energy per- 
turbations. 

In Figure[3we plot the marginalized 1 and 2a contours 
in wo — wi plane. We find that the cross-correlation 
data provide only weak upper limits on these parame- 
ters, Wo < —0.53 and wi < 0.84. If we limit our analysis 
to models with c\ = the constraints are wo < —0.82 
and W\ < 0.84 at 68%. These limits are in agreement 
with similar constraints from SN-Ia as derived in previ- 
ous analysis (see for instance |4^). 

V. DISCUSSION 

The cross-correlation between CMB and LSS data pro- 
vides a new complementary way of testing dark en- 
ergy. By isolating the ISW contribution to the CMB 
anisotropics, the correlation is a sensitive probe of the 
dark energy properties. In this paper we have studied its 
dependence on the dark energy equation of state w and 
sound speed eg. In particular we have shown that the 
redshift dependence of the cross-correlation signal may 
discriminate between different dark energy models and 
provide constraints on the dark energy parameters alter- 
native to those inferred from cosmological distance mea- 
surements. In addition the ISW correlation is insensitive 
to a number of cosmological parameters which usually 
limits the dark energy parameter estimation from CMB 
data alone. A precise knowledge of the galaxy bias is nec- 
essary for this type of measurements to be competitive. 
In fact a redshift or scale dependent bias can in principle 
mimic the effect induced by dark energy. 

We have also briefly reviewed the current observational 
status and inferred constraints on w and c\ using the cur- 
rent limited datasets. We found that, even under several 
theoretical and optimistically experimental assumptions, 
the actual constraints are weak. However, the presence 
of a dark energy component is clearly significant in the 
data and interesting. The constraints on the equation of 
state strongly depends on the dark energy sound speed. 
In the case c\ = \ the ISW data provide only a weak up- 
per limit on w < —0.53 at ler, while tighter bounds are 
obtained assuming c\ — 0. In agreement with previous 
independent works based on CMB data alone we found no 
constraints on the dark energy sound speed. Slowly vary- 
ing dark energy models are also consistent with current 
ISW data but not significantly preferred. The upcoming 
deep redshift surveys such as LSST, KAOS or ALPACA 
are optimal datasets for studying the redshift evolution 
of the CMB-LSS correlation and provide alternative con- 
straints on the dark energy parameters. 
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